Observational constraints on dark energy with generalized equations of state 
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We investigate the effects of viscosity terms depending on the Hubble parameter and its derivatives 
in the dark energy equation of state. Such terms are possible if dark energy is a fictitious fluid 
originating from corrections to the Einstein general relativity as is the case for some braneworld 
inspired models or fourth order gravity. We consider two classes of models whose singularities in 
the early and late time universe have been studied by testing the models against the dimensionless 
coordinate distance to Type la Supernovae and radio-galaxies also including priors on the shift and 
the acoustic peak parameters. It turns out that both models are able to explain the observed cosmic 
speed up without the need of phantom (w < — 1) dark energy. 



PACS numbers: 98.80.-k, 98.80.Es, 97.60.Bw, 98.70.Vc 

I. INTRODUCTION 

According to the new cosmological picture emerging 
from data (only few years ago unexpected), we live in 
a spatially flat universe filled with a subcritical matter 
content and undergoing an accelerated expansion. The 
Hubble diagram of Type la Supernovae (hereafter SNela) 
has been the first cornerstone, but quite soon other ob- 
servational data, from the cosmic microwave background 
(hereafter CMBR) anisotropy spectrum 2] to the large 
scale structure properties Q , further corroborate the first 
impression from SNela. The astonishingly precision of 
the CMBR spectrum measured by the WMAP satellite 
PI and the Gold SNela sample of Riess et al. Q repre- 
sent the last and still more convincing evidences in favor 
of this new description of the universe. 

Both the cosmic speed up and the flatness of the uni- 
verse have posed serious problems to the cosmologist 
community. Since matter alone (both visible and dark) 
cannot be enough to close the universe, a new component 
was invoked as a dominant term. Moreover, in order to 
explain the observed cosmic speed up, this new fluid must 
have a negative pressure. Being obscure both in its ori- 
gin and its properties, this component was baptized dark 
energy. Understanding its nature and nurture is one of 
the most fascinating and debated challenges of modern 
cosmology. Although the old cosmological constant 
may play naturally the role of dark energy and also fits 
well the observed dataset 0, j it is affected by serious 
theoretical shortcomings that have motivated the search 
for alternative candidates. It is nevertheless worth ob- 
serving that, to a large extent, the different models may 
be broadly classified in three classes which we will refer 
to as scalar fields, unified dark energy (hereafter UDE) 
and modified gravity. (Sometimes, some of these models 
maybe rewritten as models from another class). 

For models belonging to the first class, dark energy 



is an effective fluid originating from an ultralight scalar 
field 4>, dubbed quintessence, evolving under the action 
of a suitably chosen self - interaction potential V{(f). The 
choice of V{(j)) is crucial and different functional expres- 
sions have been investigated, from power - law 9] to expo- 
nentials 0, ^1 and a combination of both (l^j ■ Without 
entering into details (for which the reader is referred to 
^3)) only stress that such models are able to cor- 
rectly reproduce the observed data, but are affected by 
fine tuning problems and the open issue of understanding 
where the quintessence scalar field comes from. 

Notwithstanding the strong efforts made up to now, it 
is not known what is the fundamental nature of the dark 
energy so that it is also worth investigating the possibility 
that dark energy and dark matter (the other unknown 
component of the universe) are two different aspects of 
the same fluid. This is the underlying idea of models 
belonging to the UDE class. In such an approach, a single 
fluid with an exotic equation of state plays the role of 
dark matter at high densities and dark energy at low 
densities. Typical examples are the Chaplygin gas [T3 |. 
the tachyonic fluid and the Hobbit 0| model^ . 

In both classes described above, dark energy is ex- 
plained in terms of a new fluid in the framework of stan- 
dard general relativity. However, the Einstein theory of 
gravity has been experimentally tested only on scales up 
to the Solar System and hence it is far from be verified 
that its validity holds also on cosmological scales. Moti- 
vated by these considerations, different models have been 
proposed where the cosmic speed up is explained in terms 
of a matter only universe regulated by dynamical equa- 



Similar t o th e UDE models is the phenomenological inflessence 
scenario llTl where a single fluid is used to explain both inflation 
and dark energy. 
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tions that are different from the usual Friedmann ones as 
a consequence of a generalized gravity theory. Examples 
are the braneworld inspired Dvali - Gabadadze - Porrati 
model |0| and fourth order theories of gravity, both in 
the metric [H IH i22j and Palatini 24^ 15j for- 
mulation. Although such models are able to give rise 
to accelerated expansion without the need of any dark 
energy, a strong debate is still open about their compati- 
bility with the standard tests of gravity in the low energy 
limit. 



From the short (and far from exhaustive) overview pre- 
sented above, it is clear that there is much confusion. 
However, it is possible to put some order in this some- 
what chaotic situation by considering a particular feature 
of the dark energy, namely its equation of state (here- 
after EoS), i.e. the relation between the pressure and 
the energy density. Indeed, whatever is the model con- 
sidered, it is always possible to introduce a sort of dark 
energy fluid whose energy density and pressure are deter- 
mined by the characteristics of the given theory. While 
for scalar fields and UDE models, the EoS is a function 
of the energy density only, in the case of modified grav- 
ity theories, the effective EoS depends also on geometry, 
e.g. on the Hubble parameter and/or its derivatives. It is 
therefore tempting to investigate the properties of cosmo- 
logical models starting from the EoS directly and testing 
whether a given EoS is able to give rise to cosmological 
models reproducing the available dataset. A first step 
forward in this direction has been undertaken in a recent 
paper where the singularities of models assigned by 
generalized EoS have been investigated. Since the main 
interest was there in asymptotic behaviors (i.e., in the far 
past and far futu re) . any matter term was not included 
in the analysis in ^(f\ which makes it possible to analyti- 
cally solve the dynamical equations. In order to elaborate 
further on this original idea, we have to include the con- 
tribution of dust (dark/baryonic) matter. The present 
work is thus complementary to |26j since i.) we consider 
the present rather than the asymptotic states of the mod- 
els, and a.) we investigate the viability of the models and 
constrain their parameters using the observational data. 



The plan of the paper is as follows. In Sect. II, we in- 
troduce the models that we investigate, chosen for their 
interesting asymptotic behaviors and give a general the- 
oretical discussion on generalized EoS. The observational 
data, the method used to constrain the model parameters 
and the results of the analysis are presented in Sect. HI. 
Some hints on the high redshift behaviour of the models 
are presented in Sect. IV, while, in Sect.V, we qualita- 
tively discuss the issue of structure formation. We sum- 
marize and conclude in Sect. VI, while further details on 
the theoretical foundation of the generalized EoS we con- 
sider are given in two appendixes. 



II. THE MODELS 

For a spatially flat homogenous and isotropic Robert- 
son-Walker universe, the Einstein equations reduce to 
the usual Friedmann equations : 

H' = YP' (1) 

2H = ~k\p + p) (2) 

where H = a/a is the Hubble parameter, a the scale fac- 
tor and an overdot denotes the derivative with respect 
to cosmic time t. On the right hand side of Eqs.l^l and 
0, we have set = SttG, while p and p are the total 
energy density and pressure respectively. Since the ra- 
diation term is nowadays negligible, we assume that the 
universe is filled only with dust matter and a dark energy 
fluid and use subscripts AI and X to denote quantities 
referring to the first and second component respectively. 
Provided the EoS pi = Pi{pi) is somewhat given, the evo- 
lution of the energy density of the i - th fluid may be (at 
least, in principle) determined by solving the continuity 
equation : 

p^ + m{p,+p^) = Q . (3) 

This may be conveniently rewritten in terms of the red- 
shift z = 1/a — 1 (having assumed a = 1 at the present 
day) as: 

d(H ^ ?,{\ + Wi)pi 

dz 1 + z ^ ' 

with Wi = Pi/ Pi the EoS parameter^. For dust mat- 
ter, wm = 0, Eq.Q) is easily integrated to give pM — 
^MPcriti^ + z)'^ with = Pm/PctH the matter den- 
sity parameter, pcrit = "^Hq/k^, the present day critical 
density and hereon quantities with a subscript are eval- 
uated at z = 0. 

Concerning the dark energy, almost nothing is known 
about its EoS given the ignorance of its nature and its 
fundamental properties. The simplest choice 

Px =wxpx , (5) 

with wx a constant, has the virtue of leading to an inte- 
grable continuity equation thus yelding : 

PX ^nxPcr^t{l + zf^^+^^'> (6) 



^ In the following, as EoS we will refer both to the relation pi = 
Pi{pi) and to Wi indifferently. 
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with r^x the dark energy density parameter. Note that, 
because of the flatness condition, it is Ojvf + i^x = 1- 

Although a constant EoS dark energy model (also re- 
ferred to as quiessence or QCDM) nicely fits the avail- 
able data, there are two serious shortcomings in this ap- 
proach. First, there are no theoretical models predicting 
a rigorously constant EoS so that this assumption lacks 
a whatever background motivation. Moreover, fitting a 
large set of data points towards wx < — 1 as best fit (see, 
e-g-, 113) so that px +Px < 0. Models having wx < —1 
are collectively referred to as phantom models |28l | and 
have some disturbing features such as violating the weak 
energy condition px +px > and leading to a divergence 
of the scale factor a{t) in a finite time (referred to as Big 
Rip). 

It is worth noting that allowing wx to evolve with 
z does not seem to solve the problem. On the con- 
trary, fitting parametrized model independent EoS to 
the available dataset, as wx = wo + wiz [29l | or wx = 
Wo + Waz/{1 + z)P with p = 1,2 j33, still points at 
Wo < —1- Moreover, the results seem to suggest that 
Wx crosses the phantom divide, i.e. the EoS changes 
from Wx > — 1 in the past to wx < — 1 today. Recently, 
this feature has been recovered in exact models where 
inflation is matched to late-time acceleration by scalar 
phantom -non -phantom transitions [sj- Other realistic 
models are achieved by introducing also dark matter into 
the game js^- 

Here, we investigate the possibility that the dark en- 
ergy EoS depends not only on the energy density px , but 
also on the Hubble parameter H and/or its derivatives. 
We will refer to such a model as generalized EoS. As a 
preliminary remark, let us note that these same models 
were also referred to as inhomogeneous EoS in |2(j | . How- 
ever, it is worth stressing that the term inhomogeneous 
does not mean that we consider an inhomogeneous uni- 
verse (so that the RW metric may still be retained) nor 
that we consider the possibility that dark energy may 
cluster in the nonlinear regime of perturbations (as, e.g., 
in [S^). We use the term inhomogeneous since assuming 
Wx = Wx (px , H, H) introduces a viscosity term in the 
continuity equation similar to what happens in fluidody- 
namics for an inhomogeneous fluid. See Appendix A for 
an explicit derivation, while specific examples of EoS due 
to time - dependent bulk viscosity are given in, e.g., [s^ . 

In the followin g, w e will describe two classes of mod- 
els introduced in jig chosen because of their interesting 
asymptotic properties. Actually, we slightly generalize 
them and find out some degeneracies from the point of 
view of observations. Finally, we discuss the general theo- 
retical features that a given generalized EoS should fulfill 
in order to match with any observational cosmology. 



Px = Wfpx + whH'^ 



(7) 



with (wfjWn) two undetermined constants^. Eq.jTI) is a 
particular case of the more general class of dark energy 
models with EoS : 



Px 



= -PX - Ap"x - BH^f" 



(8) 



that reduces to Eq.Q by setting a = (3 = 1, wj = 
— (1 + A) and B = wh- The singularities (in the far 
past and in the far future when the matter term may 
be neglected) of this class of models have been investi- 
gated in detail in where it has been shown that a wide 
range of interesting possibilities may be achieved depend- 
ing on the values of (a, /3). In particular, for a = /3 = 1, 
even if < —1, the Big Rip is avoided provided that 
Weff = Wf + whu^ /'i > — 1 (for general classification of 
future DE singularities, see|33)- While the general model 
in Eq.Q is characterized by five parameters (namely, a, 
(3, A, B and 51 a/) and is therefore quite difficult to con- 
strain observationally, the case in Eq.© has the virtue 
of avoiding the Big Rip and (possibly) the need for phan- 
tom fields with the further advantage of being assigned 
by only three parameters. For these reasons, we will con- 
sider hereon only this particular realization of Eq. (jSJ . 

The EoS may be written in a more significant way 
using Eq.JQ) to express as function of px and pM- 
We thus get : 



PX = WfPx +WHH^[^Mil + + 7]{z)] 
wl^^ilM{l + Z) 



wde 



Px{z) , (9) 



so that the dark energy EoS reads : 



Wx = Wde 



'fi/nM{i + zf 

Tj{z) 



(10) 



In Eq.Q, in the first row, we have rewritten Eq.QJ as 

ffV^q = = ^Mil + zf + Tj{z) with T] = px/pcr^t, 

while, in the second row, we have defined : 



Wde = Wf + 



k^wh 



eff 



(11) 



A. Increased Matter 

^ Note that, while Wf is dimensionless, the dimensions of wh are 
Let us first consider the case : fixed in such a way that wgH'^ is a pressure term. 



4 



Ea. (|10|l nicely shows that the model we are considering 
is a simple generalization of the QCDM case to which it 
reduces when wh = 0. For wh 7^ 0, wx{z) may still be 
expressed analytically. To this aim, let us insert Ea. l|lU|) 
in the continuity equation for px and integrate it with 
the initial condition px{z = 0) = i^xPcrit- Using the 
dimcnsionless quantity rj ~ px/ Pcrit , we find : 



r,^(nx + ^^^m!^] (i+,)3(i+^..)„!!^m^(1+,)3 



Wde 



Wde 



Inserting now Eq. (|12|l into Eq. , we obtain : 



E'(z) = nM{i + zf + (1 - nM){i + z)3(i+---) (13) 

with 



M 



M 



eff ' 
WDE 



(14) 



We thus get a quite interesting result. The EoS (|10|) leads 
to the Hubble parameter that has formally the same ex- 
pression as that of QCDM models, but at the price of 
shifting both the matter density parameter and the dark 
energy EoS. As such, a fitting analysis based on observ- 
ables depending only on E(z) as is the case of the SNela 
Hubble diagram (see later) will give biased results. In 
particular, if wde < and wh > 0, then flM > 
so that the best fit values {■w{^*^,ft'^/) may erroneously 
point towards a universe with a higher matter content 
and a phantom dark energy. Because of this peculiar- 
ity, we will refer in the following to this scenario as the 
increased matter (hereafter IM) model, even if the in- 
creasing of the matter content is only apparent. 

It is interesting to note that, fitting only the SNela 
Gold dataset without any prior on flM, Jassal et al. 
have found SIa/ — 0.47 as best fit, significantly higher 
than the fiducial value Qm — 0.3 suggested by cluster 
abundances. On the other hand, fitting the WMAP data 
only gives ^Im — 0.32 as best fit value thus rising the 
problem of reconciling these two different estimates. In 
their conclusions, Jassal et al. argue in favor of possi- 
ble systematic errors in one of the dataset. Eas. (|13|) and 
(|14|1 . however, furnish an alternative explanation. Sim- 
ply, since (as we will see more in detail later) the SNela 
Hubble diagram only probes E{z), fitting this dataset 
with no priors gives constraints on f2jv/ rather than ^Im- 
Since flM > ^m, Eg. 1)14(1 suggests that a a possible way 
to reconcile this worrisome discrepancy may be to resort 
to our model assuming w'^/ /wde < 0. Since it is reason- 
able to expect that wde < 0, wh > have to imposed. 

Provided a suitable method is used to determine ^Im, 
Wde and CIm, it is possible to get an estimate of w'^'^ as 



w 



;// _ wde{^m — ^m) 



M 



M 



(15) 



and then of the true barotropic factor in Eg. 110(1 as : 



Wf = Wde - wll^ 



(16) 



which shows that, assuming wh > for what said above, 
leads to Wf < Wde- As a consequence, if the effective 
barotropic factor wde is in the phantom region, the true 
one Wf will be deeper in the phantom. As such, the 
model is thus not able to evade the need for phantom 
dark energy, but it is still interesting since it makes it 



(12) possible to avoid any Big Rip if wde 



eff 

Hi 



is larger 



than the critical value {wx = — 1)- 

There is, actually, an easy generalization of the model 
able to both eliminate future singularities and phantom 
fields. To this aim, let us consider the following EoS : 



Px = Wfpx + whH'^ + WdnH 



(17) 



with WdH a new constant. Ea. ((17() is inspired by Eq.(47) 
in 26] to which it reduces by setting Wf = w, wh = 
~2>{1 ~\- w) / and WdH = — 2/k^. It is easy to show that, 
with these positions and in absence of matter, Eg. l(17() 
is an identity. When a matter term is added and the 
constants (w/, wh, Wdn) are left free, this is no more true 
and interesting consequences come out. To see this, let us 
use the Friedmann equations to get and H in terms 
of (p M, Px, Px )■ It is then easy to show that Eq. ((17(1 may 
be rewritten as : 



WdHK" \ I whk"^ WdHn'^" 

1 H — \Px = \Wf ^ — ) Px 



3 2 

/ Whu"^ WdHn"^ 



\ 3 



PM (18) 



so that the dark energy EoS is given again by Eg. l(10() 
provided Eos. ((11(1 are generalized as : 



WDE 



eff 



Wf 



k^wh 



WdHl^ 



WdHl^- 



k^wh 



WdHK 



WdHK 



(19) 

while f2jvf is still defined by Eg. 1(141) with wde and w'^f/ 
now given by Egs. ((19() . By using this generalization, 
we can still get 51 jv/ > ^m, but now the conditions 
w'^j^ /wde > and wde < does not imply wh > 0. 
As a further consequence, the true barotropic factor Wf 
is now given as : 



Wf^[l + ^^]{WDE-W';j/) (20) 



and we can get Wf > —1 provided that the condition 
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WdH > ^ 



WDE 



WDE 



eff 



(21) 



is verified. Since Eqs.Q and (|18|l lead to models that 
are fully equivalent from the dynamical point of view, we 
will consider hereon these models as a single one referred 
to as the increased matter model. 



B. Quadratic EoS 

Let us now consider a different approach to the dark 
energy EoS. As yet stated above, we may think at the 
EoS as an implicit relation such as : 

F{px,Px,H)^0 

which is not constrained to lead to a linear dependence 
of px on px- As a particularly interesting example, let 
us consider the case : 



{px +Px)^ -Cspji 1 



H 







(22) 



with Cs and Hg two positive constants. Ea. H22l) has been 
proposed in '2^ where it has been shown that the cor- 
responding cosmological model presents both Big Bang 
and Big Rip singularities. As such, this model may be 
a good candidate to explain not only the present state 
of cosmic acceleration, but also the inflationary one and 
is therefore worth to be explored. Since its EoS is given 
by a quadratic equation in wx , we will refer to it in the 
following as the quadratic EoS (hereafter QE) model. 

As a first step, it is convenient to rewrite the continuity 
equation for px as a first order differential equation for 
E{z) = H{z)/Hq. To this aim, one may differentiate 
both members of Eq.([3l with p — px + ^MPcrui^ + z)^ 
and solve with respect to dpx/dz. On the other hand, 
from Eq. , we easily get : 



il + wx)px = ±px\ Cs 1 



El 

E 



(23) 



where the sign changes for wx = — 1 (phantom divide), 
corresponding to the standard ACDM model. On the 
other hand, the sign changes also for Es = E which cor- 
responds to a finite redshift z = Zs- The evolution of the 
Hubble parameter is achieved by inserting the previous 
expressions for dpx/dz and (1 + wx)px in Eq.Q, we 
finally get : 



which may be integrated numerically with the initial con- 
dition E{z = 0) = 1 provided that values of (f^M, Es, Cs) 
are given. To this end, it is convenient to express Cs in 
terms of a more common quantity. To this aim, let us 
remember that the deceleration parameter q = —ad/ci^ is 
related to the dimensionless Hubble parameter E{z) as : 



1 + z dE{z) 
E{z) dz 



(25) 



so that {dE/dz)z=o = 1 -I- go- Evaluating {dE/dz)z=o 
from Eq. and solving with respect to Cg , we get : 



Cs = 



[3»M - 2(1 + go)]' 

9{i ~ nM)Hi - Es) 



(26) 



where no sign ambiguity is present. Hereon we will 
characterize the model by the values of the parameters 
{qo,Es,flM)- It is interesting to note that Cs vanishes 
for flM = 2(1 + qQ)/5 which takes physically acceptable 
values for —1 < go < 1- In such cases, Ea. p^ reduces 
to Px — —px, i-c. we recover the usual ACDM model 
for which it is indeed qo = —1 + 3riM/2. On the other 
hand, Cs seems to diverge for Hm = 1- Actually, this is 
not the case. Indeed, for ^Im = 1, we have a matter only 
universe for which it is qo — 1/2 thus giving Cs = and 
E'^ = ilMil + z)^ so that Ea. (|^ is identically satisfied. 
Once Ea. (|24|) has been integrated, one may use Eq.Q to 
get px{z) and hence wx{z) from Ea. (|^ . 

Some important remarks are due at this point. If we 
start from a non- phantom phase (1 -I- wx > 0) impos- 
ing z = as initial condition, for large z (far past) 
the universe remains in a non -phantom phase. This is 
the case which we will match with observational data. 
On the other hand, starting from a phantom regime 
(1 + Wx < 0), we can reach a phase where 1 + wx = 
or E = Es in a finite time, i.e. for finite z. In this 
case, E{z) decreases with z, reaches Es for a finite Zs 
and after starts to increase. The sign in Ea. H23|l hanges 
at E = Es, then Es is a minimum for E(z). This means 
that, for z < Zs, the sign of the second term of r.h.s of 
Ea. (|24|) is minus, while, for z > Zs, the sign is plus. In 
conclusion, depending on the nature of dark energy fluid 
{wx), the model can represent a non- phantom cosmol- 
ogy or the evolution from a phantom to a non - phantom 
cosmology through the crossing of the phantom divide. 
For sake of simplicity, we will match with observations 
only the non - phantom solution since we want to test the 
viability of generalized EoS but, in principle, the method 
to constrain method which we will discuss below can be 
extended also to phantom cosmology. 



dE _ snuii + z)^ 

~ 2E 

, 3[E^-nM{l + z)^] [Cs{l-Es/E)]'/\ ^^^ 
2E{l + z) ^ ' 



C. Generality of generalized EoS of the universe 

The above IM and QE models are particular examples 
of generalized EoS which can be used to flt observational 
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data. A general approach to our generalized EoS is pos- 
sible from a theoretical point of view. We will show that 
any observational cosmology may emerge from such EoS. 
By using a single fimction f(t), we now assume the fol- 
lowing EoS : 



4/' / 



p 



(27) 



Then it is straightforward to show that a solution of the 
Friedmann equations (0), and the continuity equa- 
tion ^ is given by 



(28) 

Then any cosmology given hy H ~ f{t) can be realized 
by the EoS (EZJ. 

As first example, if f{t) and therefore H is given by 



Since 



a 
a 



al3 



13(1 + a)t^ 



(35) 



if i > 0, there are two solutions of a = , one corresponds 
to late time and another corresponds to early time. The 
late time solution i = of a = is obtained by neglecting 
7. One obtains 



l/p-a-l 



(36) 



On the other hand, the early time solution could be found 
by neglecting f3, which is O (10~^), to be 



m = H 



3 {w,n + 1) \< 

the EoS has the following form : 



(29) 



(37) 



Then the universe undergoes accelerated expansion when 
< t <te and ti <t<ts. Here ts is Rip time: 



P = 
Since 



3{w,n + 1)k\ p 



(30) 



1//3 



(38) 



te may be identified with the time when the infiation 
ended. One is able to define the number of the c - foldings 
Np as 



H = 



the EoS parameter WcS defined by 



P , 2H 

Wc« = - = -I - X772 ' 



(31) 



(32) 



goes to Wm > ~1 when t ^ and goes to —2 — Wm < — 1 
at large times. The crossing WeS = —1 occurs when 
= 0, that is, t = ts/2. Note that 



a 
a 



27 {Wra + lf {ts~tft^ 



X it- 



(3w„i + l)ts 



(33) 



Hence, if Wm > —1/3, the deceleration of the universe 
turns to the acceleration at t = ta = {3wm + l)is/4- 
As second example, we now consider 



fit) - H{t) 



apt 



-Ml-fln(7+^) (7+^) 



(34) 



= In (a {te) /a (0)) 



Hence, it follows 



Np = -a In 



'l^|ln(27)\ 
.l-|ln(7)y 



(39) 



(40) 



Thus, we demonstrated that generalized EoS of the uni- 
verse, in the same way as scalar model of Ref.'sT,'^ may 
present the natural unification of the early time infiation 
and late time acceleration. Working in the same direc- 
tion, with more observational constraints, one can sug- 
gest even more realistic generalized EoS of the universe, 
describing the cosmological evolution in great detail. 

However, some comments are necessary at this point. 
As for the specific models analyzed above, to test mod- 
els against observations is convenient to translate all in 
terms of z, i.e. 



H{z) = f{z,z) = - 



z+l 



(41) 



and then E{z). This allows to select interesting ranges 
of z to compare withobservations, for example 100 < 
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z < 1000 for very far universe (essentially CMBR data), 
10 < z < 100 (structure formation), < z < 10 (present 
universe probed by standard candles, lookback time, 
etc.). Then different datasets, coming from different ob- 
servational campaigns, have to be consistently matched 
with the same cosmological solution ranging from infla- 
tion to present accelerated era. This program could be 
hard to be realized in details because of the difflculties 
to join together a reasonable patchwork of data coming 
from different epochs, but it is possible in principle. 



III. CONSTRAINING THE MODELS 

Having in mind the theoretical considerations of last 
subsection, let us now develop a method to constrain dark 
energy models with generalized EoS against observations. 
In particular, we will use data coming from SNela and 
radio-galaxies for the above IM and QE models which 
reproduce the present day cosmic acceleration by means 
of inhomogeneous corrections in the dark energy EoS. 



A. The method 

In order to constrain the EoS characterizing parame- 
ters, we maximize the following likelihood function : 



C oc exp 



(42) 



where p denotes the set of model parameters and the 
pseudo - merit function is defined as : 



N 



-I 2 



i=l 

7^(p) 



1.716 



0.062 



^(p) - 0.469 
0.017 



.(43) 



Let us discuss briefly the different terms entering Ea. (|^ . 
In the first one, we consider the dimensionless coordinate 
distance y to an object at redshift z defined as : 



dz' 



(44) 



and related to the usual luminosity distance Dl as = 
(1 + z)r{z). Daly & Djorgovki |33 have compiled a sam- 
ple cornprising data on y{z) for the 157 SNela in the Riess 
et al. Q Gold dataset and 20 radio-galaxies from [ssj . 
summarized in Tables 1 and 2 of |37l | . As a preliminary 
step, they have fitted the linear Hubble law to a large set 
of low redshift {z < 0.1) SNela thus obtaining: 



h = 0.664±0.008 . 



We thus set h = 0.664 in order to be consistent with 
their work, but we have checked that varying h in the 
68% CL quoted above does not alter the main results. 
Furthermore, the value we are using is consistent also 
with Ho = 72±8 km s'^ Mpc'^ given by the HST Key 
project based on the local distance ladder and with 
the estimates coming from the time delay in multiply 
imaged quasars 40] and the Sunyaev- Zel'dovich effect 
in X - ray emitting clusters [ij • 

The second term in Ea. H43|l makes it possible to extend 
the redshift range over which y(z) is probed resorting to 
the distance to the last scattering surface. Actually, what 
can be determined from the CMBR anisotropy spectrum 
is the so called shift parameter defined as [42, H^l : 



(45) 

where zis is the redshift of the last scattering surface 
which can be approximated as |44|.: 



zis = 1048 (l + 0.00124CJ, 



-0.738 



)il + 9iLuf,) (46) 



with LOi = flih^ (with i — b, M for baryons and total 
matter respectively) and (51,32) given in Ref. 01- The 
parameter u)b is well constrained by the baryogenesis cal- 
culations contrasted to the observed abundances of pri- 
mordial elements. Using this method, Kirkman et al. 
[45I I have determined : 

LJb = 0.0214±0.0020 . 

Neglecting the small error, we thus set ojb = 0.0214 and 
use this value to determine z/g. It is worth noting, how- 
ever, that the exact value of zis has a negligible impact 
on the results and setting zis = 1100 does not change 
constraints on the other model parameters. 

Finally, the third term in the definition of takes into 
account the recent measurements of the acoustic peak in 
the large scale correlation function at 100 Mpc sep- 
aration detected by Eisenstein et al. using a sample 
of 46748 luminous red galaxies (LRG) selected from the 
SDSS Main Sample 47]. Actually, rather than the posi- 
tion of acoustic peak itself, a closely related quantity is 
better constrained from these data defined as |4^ : 



A = 



Zlrg 



Zlrg 
E{zlrg 



-i/{zlrg] 



1/3 



(47) 



with zlb,g = 0.35 the effective redshift of the LRG sam- 
ple. As it is clear, the A parameter depends not only on 
the dimensionless coordinate distance (and thus on the 
integrated expansion rate), but also on Q.m and E{z) ex- 
plicitly which removes some of the degeneracies intrinsic 
in distance fitting methods. Therefore, it is particularly 
interesting to include ^ as a further constraint on the 
model parameters using its measured value [46| : 
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Par 


bf 


la CR 


2(7 CR 


wde 


-1.03 


(-1.18,-0.91) 


(-1.39,-0.82) 




0.31 


(0.27,0.36) 


(0.23,0.41) 


Qui 


0.29 


(0.27,0.31) 


(0.25,0.33) 


iff 

■wx{z = 0) 
Qo 

ZT 

to (Gyr) 


0.018 
-1.13 
-0.54 
0.61 
14.08 


(-0.022,0.081) 
(-1.39,-0.92) 
(-0.71,-0.41) 

(0.51,0.71) 
(13.45,14.71) 


(-0.068,0.142) 
(-1.71,-0.76) 
(-0.94, -0.31) 

(0.41,0.83) 
(12.90, 15.44) 



TABLE L Summary of the results of the likelihood analysis 
for the IM model. The maximum likelihood value (bf) and 
the 1 and 2a confidence range (CR) for the model parameters 
{woe,Qm,Qm) and some derived quantities are reported. 



1.75 



1 . 5 




0.5 1 1.5 2 



FIG. 1: The best fit curve superimposed on the observed 
dimensionless coordinate distance y{z) for IM model. 



A = 0.469±0.017 . 

Note that, although similar to the usual introduced 
in statistics, the reduced (i-e., the ratio between the 
and the number of degrees of freedom) is not forced 
to be 1 for the best fit model because of the presence of 
the priors on TZ and A and since the uncertainties Ui are 
not Gaussian distributed, but take care of both statistical 
errors and systematic uncertainties. With the definition 
(|42|l of the likelihood function, the best fit model parame- 
ters are those that maximize C{p). However, to constrain 
a given parameter pi , one resorts to the marginalized like- 
lihood function defined as : 



^P^{Pi) « J dpi... J dpi^i J dpi+i... J dpnC{p) (48) 

that is normalized at unity at maximum. Denoting with 
Xo the value of the for the best fit model, the 1 and 2a 
confidence regions are determined by imposing Ax^ — 
X^ — Xo — 1 ''^^'^ = 4 respectively. 

B. Results for the IM model 

Before discussing the results of the likelihood analy- 
sis, it is worth stressing that the method described above 
makes it possible to obtain constraints on ^Im directly 
and not only on Actually, since we have used, 

in our analysis, only the dimensionless coordinate dis- 
tance, the likelihood function C should depend only on 
the parameters entering E(z) which determines the di- 
mensionless coordinate distance through Ea. H44|l . As a 
consequence, only ^Im and wde could be constrained 
without any possibility to infer Om- The use of the pri- 
ors on the shift TZ and acoustic peak parameters A is 
the way to break the degeneracy between V,m and the 
ratio w'jf/ /wjjE because both TZ and A depends explic- 
itly on ^Im and not implicitly through the (integrated or 



not) Hubble parameter. This example stresses the need 
for going beyond the SNela Hubble diagram in order to 
not only lessen the impact of eventual systematic errors, 
but also directly break possible theoretical degeneracies 
among model parameters. Table I reports a summary of 
the results giving best fit values at 1 and 2a confidence 
ranges for some interesting derived quantities. Since the 
uncertainties on {wdEi^m ,^m) are not Gaussian dis- 
tributed, we do not apply a naive propagation of errors 
to determine the constraints on each derived quantity. 
We thus estimate the 1 and 2a confidence ranges on the 
derived quantities by randomly generating 20000 points 
in the space of model parameters using the marginalized 
likelihood functions of each parameter and then deriv- 
ing the likelihood function of the corresponding quantity. 
This procedure gives a conservative estimate of the un- 
certainties which is enough for our aims. Note that this 
procedure may also be used if an analytical expression is 
not available. 

The best fit model parameters gives a theoretical y{z) 
diagram which is in very good agreement with the ob- 
served data as convincingly shown in Fig.^ It is worth 
noting that the effective dark energy EoS wde quite 
close to the ACDM value which is not surprising since 
this latter model is known to fit very well a larger set 
of observations (including also CMBR anisotropy spec- 
trum and matter power spectrum 0, H])- For the same 
reason, it is also not unexpected that flM turns out to 
be perfectly consistent with estimates coming frorn, e.g., 
WMAP CMBR spectrum ^ and clustering data %. As 
a result, the best fit (Im is only slightly larger than Hm 
thus suggesting that w'^// is vanishingly small. As re- 
ported in Table I, the best fit w'^// is indeed very small 
and w'^// = is well within the la confidence range. 
Based on Ea. lfTT)l . one may argue that wh is practically 
null and hence conclude that there is no need of any 
correction term in the dark energy EoS. However, such a 
result is biased by a theoretical prejudice. Actually, if we 
assume that Ea. (|10|l is correct, then the present day EoS 
turns out to be in the phantom region wx < — 1 and. 
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Par 


bf 


1(7 CR 


2a CR 


qo 


-0.54 


(-0.66, -0.48) 


(-0.78,-0.42) 


Es 


0.941 


(0.937,0.945) 


(0.924,0.988) 


Q,M 


0.28 


(0.26,0.29) 


(0.24,0.31) 


wx{z = 0) 


-0.93 


(-1.0,-0.71) 


(-1.0, -0.50) 


Zt 


0.70 


(0.62,0.76) 


(0.48,0.82) 


to (Gyr) 


14.21 


(13.89, 14.56) 


(13.37, 14.85) 



TABLE II: Summary of the results of the likelihood analysis 
for the QE model. The maximum likelihood value (bf) and 
the 1 and 2a confidence range (CR) for the model parameters 
(qo, Esj^Im) and some derived quantities are reported. 



using Ea. H16|l . the true barotropic factor Wf is estimated 
(at 2a CR) as : 

-1.35 <Wf < -0.75 

which is still more in the domain of phantom fields. A 
possible way to recover Wf > —1 without altering none of 
the fit results is to resort to Ea. (|18|l rather than Eq.(|5Jl. 
In such an approach, while the present day value of the 
dark energy EoS wx{z — 0) is still given by the value 
reported in Table I, the constraints above on Wf are no 
more valid and should be replaced by : 



-1.35 <Wf[l 



WdHK 



,2 



< -0.75 



It is then possible to fit the same dataset with the same 
accuracy by choosing a whatever value of Wf > —1 pro- 
vided that wj^H is then set in such a way that the above 
constraint is not violated. 

Having discussed the constraints on the dark energy 
EoS, it is interesting to compare some model predictions 
with the current estimates in literature. As a first test, 
we consider the present day value qo of the decelera- 
tion parameter. Since, as a general rule, go = 1/2 -f 
{3/2)flxWx{z — 0) and we get that both = 1 ^ 
and Wx{z = 0) are consistent with the ACDM values, it 
is not surprising that the best fit qo = —0.54 is quite 
similar to the concordance model prediction go — —0.55 
^lEQiHl' On the other hand, go may also be determined 
in a model independent way by using a cosmographic ap- 
proach, i.e. expanding the scale factor a{t) in a Taylor 
series and fitting to distance related data. Using a fifth 
order expansion and the Gold SNela sample, John 
has obtained go = — 0.90±0.65 which is a so large range to 
be virtually in agreement with almost everything. How- 
ever, we note that our best fit value is significantly larger 
than his estimate. Although this could be a problem, we 
do not consider it particularly worrisome given the uncer- 
tainties related to the impact of truncating the expansion 
of a{t) to a finite order. 

According to the most recent estimates, the cosmic 
speed up has started quite recently so that it is possi- 
ble to detect the first signature of the past decelerated 
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FIG. 2: The best fit curve superimposed on the observed 
dimensionless coordinate distance y{z) for the QE model. 



expansion. This opens the way to determine the transi- 
tion redshift defined as the solution of the equation 
q{zT) — 0. For the lAI model, it is possible to derive an 
analytical expression : 



ZT = 



1 - 2go + 3wde 
(1 - 2qo){l + WDE) 



1 . 



(49) 



The constraints we get are summarized in Table I. As a 
general remark, we note that the best fit value zt = 0.61 
is once again close to that for the concordance ACDM 
model for the same reasons explained above for go. Most 
of the zt estimates available in literature are model de- 
pendent since they have been obtained as a byproduct 
of fitting an assumed model (typically the ACDM one) 
to the data. A remarkable exception is the reported 
zt = 0.46±0.13 found by Riess et al. Q fitting the ansatz 
q{z) = qo + {dq/ dz)z=oz to the Gold SNela sample. Al- 
though our best fit estimate is formally excluded at la 
level by the result of Riess et al., there is a wide overlap 
between our estimated la CR and that of Riess et al. so 
that we may safely conclude that the two results agree 
each other. 

As a final check, we consider the present age of the 
universe that may be evaluated as : 



to = tu 



dz 



{^ + ^)E{z) 



(50) 



with tH = 1/Hq = 9.778/1-1 q^^. ^-^^ Hubble time. Our 
best fit estimate turns out to be to — 14.08 Gyr with the 
2a CR extending from 12.90 up to 15.44 Gyr. This result 
is in satisfactory agreement with previous model depen- 
dent estimates such as to = 13.24^q4^ Gyr from Tegmark 
et al. and to = 13.6±0.19 Gyr given by Seljak et al. 
|j8|. Aging of globular clusters [43| and nucleochronol- 
ogy |50j give model independent (but affected by larger 
errors) estimates of to still in good agreement with our 
one. Note also that, the larger value of the best fit to we 
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FIG. 3: The percentage difference /S.H with respect to the 
ACDM scenario for the IM model. We set the parameters 
9.M (for both the IM and the ACDM model) and Q.m to 
their best fit values in Table I, and choose three values for 
wde, namely wde = —1.39 (short dashed), —1.03 (solid) 
and —0.82 (long dashed). 



FIG. 4: The percentage difference /S.H with respect to the 
ACDM scenario for the QE model. We set the parameters 
Q.M (for both the QE and the ACDM model) and go to their 
best fit values in Table II, and choose three values for Es, 
namely Es = 0.924 (short dashed), 0.941 (solid) and 0.988 
(long dashed). 



obtain is also a consequence of our assumption h — 0.664 
that is smaller than h — 0.71 used by Tegmark et al. 
and Seljak et al. as a result of their more elaborated and 
comprehensive data analysis. 



C. Results for QE model 

The results of the likelihood analysis obtained in this 
case are summarized in Table II where constraints on 
the model parameters (go, Es,Qm) and on some derived 
quantities are given, while Fig.[21 shows the best fit curve 
superimposed to the data. As it is clear, the agreement 
is quite good and, to a large extent, the QE model works 
as well as the IM one in fitting the used dataset. 

Although the dark energy EoS is significantly differ- 
ent, it is interesting to note that most of the QE model 
predicted quantities are in very good agreement with the 
same quantities for the IM scenario. It is worthwhile to 
spend some words on the present day value of the deceler- 
ation parameter. While for the lAI model qo is a derived 
quantity, here qo is used to assign the model characteris- 
tics and is therefore obtained directly from the likelihood 
procedure. Notwithstanding the different role played in 
the two cases, the best fit value is identical, while the 
1 and 2a CR overlap very well'*. This result may be 
qualitatively explained noting that, in both cases, set- 
ting the model parameters to their best fit values leads 
to a cosmological model that is quite close (from a dy- 



* It is not surprising that the 1 and 2a CR on qo in Table I are 
slightly larger than those in Table II. Actually, for the IM model, 
the constraints on qo are determined by the uncertainties on 
the other model parameters, while, in the present case, they are 
obtained from the fitting procedure directly. 



namical point of view) to the concordance ACDM sce- 
nario. As a further support to this picture, let us note 
that wxiz = 0) ~ — 1 as for the cosmological constant. 
Note that, because of our choice to be in a non-phantom 
regime, wx is forced to be larger than —1 so that we may 
give only upper limits on this quantity. If we had cho- 
sen the phantom regime, and the possibility to cross the 
phantom divide, we should have given also lower limits 
to Wx- 

The fact that the model is indeed close to the usual 
ACDM can be also understood looking at the constraints 
on Eg. Eg. lf^ shows that the ACDM model is obtained 
in the limit E^ = 1 so that, in a certain sense, this pa- 
rameter measures how far the dark energy fluid is from a 
simple cosmological constant. Indeed, the best fit value 
i?s = 0.941 is quite close to 1 so that one could argue that 
only minor deviations from the standard ACDM model 
are detected. 

As a final remark, we note that both the transition red- 
shift zt (that, in this case, has to be evaluated numer- 
ically) and the age of the universe to are in agreement 
with what is expected for the ACDM model (once the 
value of h is raised to the commonly adopted h = 0.70) 
so that it is not surprising that we find good agreement 
with the values reported in literature and quoted above 
when discussing the results for the IM model. 

IV. HIGH REDSHIFT BEHAVIOUR 

The likelihood analysis presented above has demon- 
strated that both the lAI and the QE model are able 
to fit the data on the dimensionless coordinate distance 
to SNelA and radio -galaxies. As a further test, we have 
also imposed priors on the shift and acoustic peak pa- 
rameters in order to better constrain the models. Nev- 
ertheless, the bulk of the data we have considered only 
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probes the low redshift {z < 1.7) epoch of the universe 
expansion. It is therefore worth investigating what are 
the consequences of the bulk viscosity corrections to the 
dark energy EoS at higher redshifts. We can infer some 
interesting qualitative results by simply comparing the 
high redshift expansion rate to those predicted by the 
successful concordance ACDM model. 

To this aim, we show in Figs.|31and0]the percentage 
deviation AH = \{H\ — H,nod) / Ha\ of the Hubble rate 
H{a) predicted by the IM and QE models respectively 
with respect to the ACDM scenario. As it is apparent, in 
both cases, AH significantly differs from the null value 
only in the very recent epochs, i.e. on the redshift range 
probed by the SNelA data. 

Actually, these results may be qualitatively explained 
for both models. As Ea. (|13|l shows, the Hubble param- 
eter for the IM model is formally the same as that of 
the ACDM model provided the matter content ^Im is 
shifted upwards being replaced by ^m- According to the 
constraints discussed above, however, and Q,m are 
almost the same and, furthermore, wde is only slightly 
different from the cosmological constant value. It is thus 
not surprising at all that the IM and the ACDM models 
predict almost the same expansion history. A similar dis- 
cussion also holds for the QE case. As yet said, the QE 
model reduces to the ACDM one in the limit Eg ^ 1 and 
indeed, the la CR for this parameter is not too far from 
1 to give rise to significant deviations at high redshift. 

Although encouraging, both these results are not sur- 
prising. Indeed, any reasonable dark energy model 
should predict that the energy density of this component 
becomes vanishingly small in the far past in order to al- 
low the universe undergoing a matter dominated phase. 
Moreover, the dark energy contribute to the total en- 
ergy budget must be negligible during the radiation dom- 
inated period so that the nucleosynthesis can take place 
in the usual way. As such. Figs . El and ^ are reassuring 
since they show that both the IM and QE models be- 
have as the ACDM one during the matter and radiation 
dominated epoch thus making us confident that the con- 
straints from nucleosynthesis are satisfied. 



V. STRUCTURE FORMATION 

Including viscosity terms in the dark energy EoS could 
impact also the formation of cosmological structures. It 
is common to investigate the evolution of matter density 
perturbations S — Spm/pm, in the linear regime, solving: 



S + 2HS-4:nGpMS^0 . 



(51) 



This relation assumes that the dark energy does not clus- 
ter on small scales so that px concurs to determining 
the evolution of S only through the Hubble parameter. 
Moreover, the derivation of Ea. (|51|) implicitly assumes 
that the General Relativity in its Einsteinian formula- 



tion is the correct theory of gravity (see, e.g., (Sll for a 
textbook description). 

The analysis of our generalized EoS presented up to 
now mainly rely on the dynamics of the cosmic expansion 
which can be investigated also assuming that both the 
lAI and QE models are phenomenological descriptions of 
an unspecified theory. Put in other words, we do not need 
to know what is the Lagrangian giving rise to Eqs.(^ and 
(0) with the dark energy EoS given by Eqs.flZI) and 
In order to investigate the issue of structure formation, 
it is, on the contrary, mandatory to have a Lagrangian 
description of both models. 

Let us first consider the IM model in which case the 
EoS (in its more general form) is given by Eq. ljlTI) . In 
Appendix B, we show that such a model could be derived 
by the following action : 



S = d Xy/^< 



2^2 



Such an action is similar to that of most of scalar field 
quintessence models with V{(j)) the self- interaction po- 
tential, but now the sign of the kinetic energy is not fixed 
from the beginning. It is possible to show that this choice 
allows to cross the phantom divide line in a natural way 
which is not the case in the usual quintessence models. 
For what we are concerned here, it is worth stressing 
that such an action makes it possible to interpret the 
IM model in the framework of Einsteinian General Rel- 
ativity"". In this leading role in determining the 
clustering properties of the dark energy is played by the 
the EoS and sound speed Cg. Assuming the pressure 
perturbations are adiabatic, the sound speed reads : 



'X 



dp 



wx + \ --r 



x 



1 d^ 

rj dz 



^ dw} 



dz 



(52) 



with rj — px{z)/ Pcrit- For the IM model, after some 
algebra, we get : 



[(1 + z)-^T] - (»M - »Af)](l + Wde)wde 
(1 + u'i3_e)(1 + z)-^r] - (Om - ^m)wde 



(53) 



with : 



= {nM-nM)ii+zf+ii-nM){i+zf<^'+^''^^ . (54) 



M — 



0, we 



It is immediately clear that, in the limit w 
get flM = and both the EoS and the sound speed 
reduces to that of a quiessence model with constant EoS, 



For a different interpretation, see |52| wliere it is sliown that 
may be derived adding a buik viscosity term in tfie en- 
ergy-momentum tensor. 
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i.e. we get = wx = wde- Indeed, for the best fit 
parameters, departs from the constant value woE only 
marginally since w)^./ is quite small. Moreover, if we 
consider the range wde < — 1, which is compatible with 
the constraints in Table I, one may easily check that, for 
z >> 0, it is : 

i.e. the IM model reduces to a phantom one with con- 
stant EoS Wx = Wde- It is easy to show that the same 
approximate result for the sound speed also holds in the 
case Wde > —1- These analogies make us confident that, 
in the IM scenario, the evolution of matter density per- 
turbations could be studied as in standard quiessence 
models (see, e.g., 53} and references therein), but dis- 
cussing further this issue is outside our aims. 

Let us now consider the QE model which, as shown in 
Appendix B, could be derived by the same action above 
choosing a different expression for uj{(j)). As a conse- 
quence, we can still use the standard theory of perturba- 
tion. Starting from Ea. (|22|l . we get : 

2^ C,E,E'[E^ -nM{l + zf] 

""^ 2{l + wx)E^[2EE'-?,nM{l + zY] ^ ' 

with Wx and Cg given by Eas. H23() and (|26(l respectively, 
E = H{z)/Hq and the prime denoting derivative with 
respect to z. Integrating Ea. l(^ for a given values of 
the model parameters inserting the result into Ea. H55|) 
above, we find out that c^(z) has an unpleasant behaviour 
since it starts from values larger than 1, which is not 
physically reasonable®, and then takes negative values 
for z > Zc, with Zc in the range (0.5,1.5) depending on 
(qo, -Es, JIm)- Positive values of the sound speed may 
be obtained also for the generalized Chaplygin gas and 
give rise to exponential blow up of the perturbations in 
the linear regime [H^l- However, it has been suggested 
that considering entropy perturbations could solve the 
problem [E^. Even if the QE model significantly differs 
from the generalized Chaplygin gas, it is possible that a 
similar mechanism could alter the sound speed in such a 
way to eliminate the disturbing feature of superluminal 
values. Indeed, when entropy perturbations enter the 
game, pressure perturbations are no more adiabatic so 
that Ea. (|52|l does not hold anymore. Whether this can 
also regularize the sound speed for the QE model is a 
topic worth to be investigated, but it is outside our aims. 

VI. CONCLUSIONS 

Dark energy EoS with inhomogeneous terms coming 
from geometry (e.g. H, H) can yield cosmological mod- 



Remember that we are working in units with c = 1 so that 
Eg. 1521 actually refers to c^/c^. 



els capable of avoiding shortcomings coming from coin- 
cidence problem and a fine tuned sudden passage of the 
universe from the decelerated regime (dark matter dom- 
inated) to the today observed accelerated regime (dark 
energy dominated). Furthermore, such models allow to 
recover also early accelerated regimes with the meaning 
of inflationary behaviors. In this paper, we have con- 
strained two physically relevant models using SNela and 
radio-galaxies data including priors on the shift and the 
acoustic peak parameters. 

In the first class of models, referred to as IM model, 
the dark energy pressure depends linearly not only on 
the energy density, but also on and H. Analytical ex- 
pressions for px{z), wx{z) and the dimensionless Hubble 
parameter E{z) have been obtained. As an interesting re- 
sult, we obtain for E{z) the same functional expression 
as for the QCDM model provided the matter density pa- 
rameter iljv/ and the dark energy EoS wx are replaced 
by biased quantities Q.m and wde- As a consequence, 
fitting to observables only depending on E{z) may lead 
to wrong results. In the second case, dubbed QE model, 
the EoS is not given explicitly, but as the solution of 
a quadratic equation. In particular, we have fitted the 
non-phantom regime wx > — 1. 

The likelihood analysis presented in Sect. HI has con- 
vincingly shown that both models are able to fit well the 
data on the dimensionless coordinate distance to SNela 
and radio - galaxies also taking into account the priors on 
the shift and acoustic peak parameters TZ and A. More- 
over, they predict values of the deceleration parameter go; 
the transition redshift zt and the age of the universe 
that are in agreement with most of the estimates available 
in literature. It is worth noting that, setting the model 
parameters to their best fit values, both the IM and QE 
models reproduce the same dynamics as the concordance 
ACDM scenario. Nevertheless, the underlying physical 
mechanism is radically different so that both approaches 
worth be further investigated. 

On the one hand, we have here only considered the 
general dynamics of the universe, i.e. the Hubble pa- 
rameter II{z), which enters in the determination of dis- 
tance related quantities. It is therefore mandatory to 
enlarge our attention considering other observational 
tests which probes the EoS itself rather than its inte- 
grated value. Ideal candidates are the growth index and 
the evolution of structures since wx{z) enters directly 
into the corresponding equations. Moreover, the full 
CMBR anisotropy spectrum has to be evaluated which 
is likely to offer the possibility to determine separately 
{wf,WH,WdH) in the case of the IM model. Actually, 
since our EoS are inhomogeneous, viscosity - like terms 
will appear in the equations governing the evolution of 
density perturbations thus probably leading to some dis- 
tinctive signatures that could be detected. 

On the other hand, further developments on the theo- 
retical side are possible. Actually, Eas. ljlTII and (|22|l have 
been assumed rather than obtained from a fundamental 
theory. As shown in Appendix B and more deeply dis- 
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cussed in [2g , it is possible to work out a modified gravity 
theory that leads to an effective dark energy fluid with 
EoS given by Ea. lfTTjl or by Ea. f^ . It should therefore 
be interesting to investigate what constraints may be im- 
posed on deviations from Einstein general relativity by 
the likelihood analysis presented here. 

As it was demonstrated in '2^, including more gen- 
eral viscosity term in the dark energy EoS brings new 
possibilities to construct both the early and late time 
universe. In particular, one may look for models that 
are able to give rise to both an early inflationary epoch 
and a present day accelerated expansion j82| . Testing 
such models against observational data is the key to se- 
lect among the different possibilities which is the most 
realistic one. 

Finally some comments are necessary on the general 
philosophy of the classes of models studied here and on 
the constraints coming from the observations capable of 
restricting and selecting viable cosmologies. Consider- 
ing also the general discussion in Sect. II (Subsect. C), 
it is clear that generalized EoS where further terms are 
added to p — —p are preferable for the following rea- 
sons. A huge amount of observational evidences indicates 
that ACDM {p = —p) is the cosmological scenario able 
to realistically describe better the today universe. Any 
evolutionary model passing from deceleration (dark mat- 
ter dominance) to acceleration (dark energy dominance) 
should consistently reproduce, based on observations, 
such a scenario. In the IM case which we have stud- 
ied, adding terms in Hubble parameter and its derivative 
allows always a comparison with standard matter param- 
eters {wm and Om) which are, in some sense, "directly" 
observable. Then the number of arbitrary choices (for ex- 
ample, to fix priors) is not so large. On the other hand, 
QE models, or in general implicit EoS, need several arbi- 
trary choices which could result completely inconsistent 
to further and more refined observations. For example, 
in our case, we have arbitrarily discarded phantom-like 
regime which could result consistent with observations 
at large z (far distances and early universe) and imposed 
arbitrary constraints on go- Due to these reasons, from 
an observational point of view, it is preferable to study 
models which imply corrections the A EoS rather than 
giving EoS in implicit form. 



APPENDIX A: BULK VISCOSITY 



Let be the four velocity of the cosmic fluid. Then 
the scalar expansion is given hy 9 = J7^^ and we have 9 = 
3H. When there is a bulk viscosity (, which could be a 
function of the energy density px of the fluid: C = C (Px), 
the first FRW equation is not changed (3/K^)i?^ — p but 
the second FRW equation and the conservation of the 
energy are expressed as 



2H + 3H^ +px-Cd 



= ^ (2ij + 3H^) +PX- 3i?C (fx) , (Al) 

= Px + 9ipx+Px-C0) 

= PX + 3H {px +PX- 3iJC {px)) . (A2) 



Then we obtain the effective pressure p as 



Px ^Px - 3HC {p 



X) 



(A3) 



Then the second FRW equation and the conservation of 
the energy have the standard forms: expressed as 



1 

PX 



2H + 3H' 



Px , 



3H {px +px) 
Then if px and px satisfy the EOS 

F{px,Px) = , 
we obtain effectively generalized EOS: 



(A4) 
(A5) 



(A6) 



F{px,Px,H) = F {px ,px + iHC{px))^0 . ( A7) 

APPENDIX B: THEORETICAL FOUNDATION 

Here, we give two possible theoretical foundations for 
the IM and QE models discussed in the text. We may 
start with the following action: 



S 



(Bl) 

Then the energy density p and the pressure p for the 
scalar field (f) are given by 



1 



2 + y(^) , p^^u; 



y(0) . (B2) 



Since we can always redefine the scalar field <^ as (/) — > 
F{(f>) by an arbitrary function F{(j)), we can choose the 
scalar field to be a time coordinate; (j> = t. Furthermore 
we consider the case that uj{4>) and V{(t)) are given by a 
single function f{(j>) as 



Then 



c.(0) = --/'(</.) , 

V{^) = 1 (3/(0)2 + /'(^)) 



p = ^m' > 



(B3) 



(B4) 
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Since p = f ^ {^n^fpf^, we reobtain the generalized 
EoS 



As an example, the case may be considered 

/ - + /i , 

with constants /o and /i. Then (|B5|I gives 



= (P + pf-^/(l-f 



(B5) 



(B6) 



(B7) 



If we neglect the contribution from the matter, from the 
first FRW equation, we find 



-p = H^ 



One may rewrite ljB7|) as 



„ . , x2 4/0 [p ( fl 



which may have a structure similar to Eq. 
also consider the case that 



/o 



which gives 



(B8) 

(B9) 
One may 

(BIO) 
(Bll) 



Hence, we obtain homogeneous equation of state with 
w ^ —1 + 2/3/0- In case without matter, by using the 
first and second FRW equations, we can rewrite (|B11|) as 



P 



1 /> + 



2(1 ~ /3) ^2 ^ 2(1 -g) 



H . (B12) 



Here a and /3 are constants. Then by identifying 

P 1 2(1-/?) 2(1 -g) 



3/oa 



(B13) 

one gets Ea. l(T7|l . while we get Eq.© for a = 1. 

We may include matter with constant EoS parameter 
Wjn = Pm/pm to the action (|B1|) . Here p^ and are 
energy density and pressure of the matter. When the 
matter is included, the FRW equations give 



= - (p,„ + pjn) , 

V{^) = ^ (SH^ + h)- P^- P"- . (B14) 



By using the energy conservation p,n+5H (p™ + Pm) = 0, 
we find Pm = PmoO'~^^^~^^'"^ with a constant pmo- Then 
if 



= — - (''Wm + 1) goe 
y(</>) = l(35'(^)2+5"(0)) 



-3(1- 



-3(l+to„)s(</,) 



(B15) 



we find a solution 
H = g\t) 



a =^ aoC' 



,ff(t) 



ao 



50 



(B16) 



Then wc find 

P -- 
P -- 



A5'(0)2_g„e-3(i+»™bW) 



1 



K 

'WmgoG 



(3.g'(0)2 + 2.g"(0)) 

3(l+«)„)s(0) 



(B17) 



Then by deleting (p from the above equations (jB17|) . we 
obtain an EoS corresponding to (jBSp . 

We may also consider the possibility to obtain the gen- 
eralized EoS from the modified gravity. As an illustrative 
example, the following action is considered: 



S = I d X^J~^ ( + -Cmattcr 



— i? + £mattcr + /(^) 



(B18) 



Here /(i?) can be an arbitrary function of the scalar cur- 
vature R and £mattcr is the Lagrangian for the matter. 
In the FRW universe, the gravitational equations are: 



= -H^ + p- f {R = 6H +12H 



+6 [ H + - H 



^ 



dt 



xf i^R = 6H + 12H' 

(2H + 3H^^ +P + f (R = 6H + UH' 



(B19) 



+2 -H -'3H'^ 
xf (r = 6i/ + 12ij2 



d^ ^^d 
d^ + ^^Jt 



(B20) 



Here p and p are the energy density and the pressure 
coming from £matter. They may satisfy the equation of 
state like p = wp. One may now define the effective 
energy density p and p by 



p- f(R^ 6H + 12H' 



+6\H + H'-H- 
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X f (^R = 6H + UH'^^ , (B21) 

p = p + J [r = &H + UH'^'^ 

X f (^R = 6H + 12H^^ . (B22) 
Thus, it follows 

p = wp+ {l + w)f (^R = 6H +12H^^ 

{ . (p 

+2 (-1 - 3w) H -'S{1 + w) H^ + -772 



+ (2 + 3w) /' [R = m+ UH^^ (.B23) 

Especially if we consider the case that / ~ foR with a 
constant /o, we obtain 

p^wp + 6/oij + 6(1 + w)foH^ , (B24) 

which reproduces I17|) by identifying p = px, P = Px , 
w = wx, 6/o = WdH 6(1 + w) — wh- Furthermore if we 
can neglect H compared with , we also obtain Eq. Q . 

It is thus demonstrated that the generalized EoS con- 
sidered in this paper maybe derived also from scalar - 
tensor or fourth order gravity theories. 
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